A consistent interpretation of the low temperature magneto-transport in graphite 
using the Slonczewski— Weiss— McClure 3D band structure calculations 
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Magnetotransport of natural graphite and highly oriented pyrolytic graphite (HOPG) has been 
measured at mK temperatures. Quantum oscillations for both electron and hole carriers are observed 
with orbital angular momentum quantum number up to N m 90. A remarkable agreement is 
obtained when comparing the data and the predictions of the Slonczewski-Weiss-McClure tight 
binding model for massive fermions. No evidence for Dirac fermions is observed in the transport 
data which is dominated by the crossing of the Landau bands at the Fermi level, corresponding to 
dE/dkz = 0, which occurs away from the H point where Dirac fermions are expected. 

PACS numbers: 



Recently, massless Dirac fermions have been observed 
at the K point of the Brillouin zone in graphene, a hexag- 
onally arranged carbon monolayer with quite extraor- 
dinary properties [JJ. Historically, graphene forms the 
starting point for the Slonczewski, Weiss and McClure 
(SWM) band structure calculations of graphite [5] . In 
graphite, the Bernal stacked graphene layers are weakly 
coupled with the form of the in-plane dispersion depend- 
ing upon the momentum k z in the direction perpendic- 
ular to the layers. The carriers occupy a region along 
the H — K — H edge of the hexagonal Brillouin zone. 
At the K point (k z — 0), the dispersion of the electron 
pocket is parabolic (massive fermions), while at the H 
point (fc z = 0.5) the dispersion of the hole pocket is lin- 
ear (massless Dirac fermions) . A clear signature of Dirac 
fermions at the H point of graphite has recently been 
reported using far-infrared magneto- absorption measure- 
ments |4 . Such measurements probe the very close vicin- 
ity of the H and K points where there is a maximum in 
the joint density of states. 

The SWM model, which provides a remarkably ac- 
curate description of the band structure, has been ex- 
tensively tested using Shubnikov de Haas, de Haas van 
Alphen, thermopower and magneto-reflectance measure- 
ments to caliper the Fermi surface of graphite [5J El [7l 
El [10] . There are even reports of a charge density wave 
state above B = 22 T P3 [TJ1 [TJ]. However, the obser- 
vation of massless carriers with a Dirac-likc energy spec- 
trum, using magneto-transport measurements [TH [TS] re- 
mains controversial, since in the SWM model, the elec- 
trons and hole carriers at the Fermi level are both massive 
quasi-particles. 

In this Letter, we report magneto-transport mea- 
surements of natural graphite at very low temperature 
(T s» 10 mK). Due to the low temperatures used, the 
magneto-transport is much richer than previously pub- 
lished data o m h m m rm rm na rm ns]. Quantum 

oscillations are observed for both majority electrons and 
holes with orbital quantum number up to almost N=100. 
We show that these oscillations are fully consistent with 



the presence of majority electron and hole pockets within 
the three dimensional SWM band structure calculations 
for graphite. At high magnetic fields (B > 2 T), a signifi- 
cant deviation from 1 /B periodicity occurs due to the well 
documented movement of the Fermi energy as the quan- 
tum limit is approached [HE]. This seriously questions 
the validity of using the high field data to extract the 
phase of the Shubnikov de Haas oscillations, and hence 
the nature of the charge carriers [TS] . 

For the measurements mm-size pieces of natural 
graphite and highly oriented pyrolytic graphite (HOPG), 
a few hundred microns thick, where contacted in an ap- 
proximate Hall-bar configuration using silver paint. The 
measurements were performed with the sample placed di- 
rectly in the mixture of a He 3 /He 4 dilution fridge, using 
an ac current of ~ 10 /iA at 10.7 Hz and conventional 
phase sensitive detection. The magnetic field was applied 
along the c-axis of the sample. 

Typical low temperature data for R xx (f2) as a func- 
tion of the magnetic field from B = — 10 T for natural 
graphite, is shown in Fig.[TJa). R XX (B) increases roughly 
linearly with the magnetic field and at B = 10 T, it is 
about three orders of magnitude larger than the zero- 
field value [T7J [TBI EH HD] • In addition, quantum oscilla- 
tions are superimposed on the large magneto-resistance 
background. These oscillations, can be better seen in the 
background removed data AR XX plotted in Fig.[TJa-c) for 
successively slower sweeps in order to reveal the quan- 
tum oscillations in the different magnetic field regions. 
The background can either be removed by subtracting 
a smoothed (moving window average) data curve or by 
numerically calculating the second derivative d 2 R/dB 2 . 
Both techniques give similar results and here we use av- 
eraging to remove the background. As the oscillations 
are periodic in 1/B, the optimal number of points used 
in the averaging depends upon the magnetic field re- 
gion. For this reason, the amplitudes of the oscillations 
in Fig. [lja-c) should not be compared, as different aver- 
aging was used to remove the background. HOPG (not 
shown) presents almost identical oscillations, with very 
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FIG. 1: (color online) (a) Right axis: Resistance R xx versus 
B measured at T = 10 mK for natural graphite, (a-c) Left 
axis: Background removed data AR XX showing quantum os- 
cillations measured over different magnetic field regions. The 
high field electron (e) and hole (h) features are indicated. The 
vertical arrows indicate spin split electron and hole features. 



slightly different frequencies, and a significantly reduced 
amplitude [21] . For this reason we concentrate here on 
the data for natural graphite. 

In the AR XX (B) data shown in Fig. [I] two series of 
oscillations can be distinguished. The oscillations start 
at a magnetic field B ss 0.1 T, and spin splitting of 
the features (indicated by arrows) is observed for mag- 
netic fields B > 1 T, compared to previous work [7] in 
which spin-splitting was only observed for the three last 
features at magnetic fields B > 2 T. The electronic g- 
factor g s can be estimated from the magnetic fields at 
which spin splitting occurs, (B z ~ 1 T), and at which 
the Shubnikov de Haas oscillations start (B c ~ 0.07 T). 
For a Landau level broadening T, spin splitting occurs 
when g s [iBB z ss T, and Shubnikov de Haas oscillations 
occur when HeB c /m* sa T where m* = 0.056m e is the 
electron effective mass [52]. Assuming T to be field in- 
dependent we can write g s « heB c /m* psB z w 2.5. 
The mobility, estimated from the condition pB c = 1, is 
~ 140, 000 cm 2 /Vs. In the Fourier transformation of the 
0-0.4 T AR XX versus (1/B) data, shown in Fig.[2ja), two 
frequencies are found and assigned to the electron pocket 
at the K point (Bp e — 6.14 T) and hole pocket at the H 
point (Bp h = 4.51 T). This assignment is the well estab- 
lished in the literature [7] [S] and it is the only assignment 
which is consistent with the magnetoreflectance measure- 
ments jTUj. For HOPG (data not shown) we find slightly 
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FIG. 2: (color online) (a) Fourier transform of the low mag- 
netic field AR XX (1/B). (b-c) Orbital angular momentum 
quantum number N, as a function of the reciprocal magnetic 
field positions of the electron and hole features, (d) and (e) 
Contour plot of the phase shift function K(<p, B) in the vicin- 
ity of the hole and electron features. Maxima in K{ip, B) 
determines the frequency and phase of the oscillations. 



higher frequencies, Bp^ = 6.49 T and Bp h = 4.73 T. 

In graphite, we have experimentally p xx 3> p xy so 
that the tensor relation for conductivity simplifies to give 
a xx oc p~ x . Therefore, conductivity maxima which occur 
at coincidence of Landau bands and the Fermi energy 
Ep [23], correspond to minima in AR XX (B). We per- 
form a classical 1 /B analysis of our data assigning an or- 
bital quantum number N to the electron and hole minima 
of AR XX (B). For N < 25 the magnetic field positions of 
each series of oscillations can be determined directly from 
AR XX (B). For N > 25, pass band frequency domain fil- 
tering was used to separate the superimposed electron 
and hole features. The position of the features in in- 
verse magnetic field versus N is shown in Fig. [2^b) . For 
both electrons and holes we can see features with angular 
quantum number 1 < N < 90 (to almost N = 100 for 
electrons) . N versus 1/B has a linear dependence and the 
slope gives the fundamental fields Bp c — 6.10 ± 0.05 T 
and Bp h = 4.50 ± 0.05 T, in good agreement with the 
values obtained from the Fourier transform. 

At low magnetic fields, i.e. at high quantum number N, 
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a perfect linear behavior in N(l/B) is observed for both 
electrons and holes. For high magnetic fields, i.e. for low 
N, clear deviations from the linear behavior are observed 
for the electron features (see Fig. [2|c)). This deviation 
from a periodic in 1 /B behavior at high magnetic fields 
is due to the Fermi level moving as the quantum limit 
is approached in graphite [5J \W\ . Clearly, the high field 
data should not be used to extract the phase of the oscil- 
lations [15 . Equally, extrapolating the low field data to 
find the intercept does not give a reliable estimate of the 
phase. Instead, we prefer to use the phase shift analysis 
method developed by Luk'yanchuk and Kopelevich [13] , 
to extract the phase from the complex Fourier transform 
f(B) of the low magnetic field AR xx (l/B). The phase 
shift function K(ip, B) = Re\e~ %v f(B)] has maximum in 
the ip — B plane which can be used to extract both the 
frequency (B) and phase (ip) of the oscillations. K(ip,B) 
is plotted in Fig. |2^d-c) in the regions of the hole and 
electron features. From the maxima, the determined fre- 
quency and phase are Bf h = 4.51 T, ipi 7 = — (0.56±0.1)7r 
and B fe = 6.14 T, ip e = -(0.86 ± 0.1)tt for the hole and 
electron features respectively For HOPG a similar anal- 
ysis gives ip h = -(1.04±0.1)tt and ip e = -(0.92 ± 0.1)tt. 

The oscillatory conductivity can be written Act cx 
cos(2irB f / B — 2irj+5) with 7 = 1/2 for massive fermions 
and 7 = for massless Dirac fermions [23l [24] . At low 
magnetic fields inter Landau level scattering is expected 
to dominate so that 6 = ir/4 for a 3D corrugated Fermi 
surface (5 — for a 2D cylindrical Fermi surface). The 
expected value of the phase tp = — 27r7 + 6 for massive 
3D fermions (7 = 1/2) is therefore ip = — 0.757T in rea- 
sonable agreement with the experimental phase for both 
electrons and holes. For HOPG the phase is also consis- 
tent with 7=1/2 but with 5 w 0. The value of 7 = 1/2 
is in agreement with published results [5] [9] and theoreti- 
cal considerations [21] ■ In contrast, the prediction for 2D 
massless Dirac fermions (7 = 0) with ip = is completely 
inconsistent with the determined phase for both electrons 
and holes. We therefore conclude that there is no evi- 
dence from transport measurements for the existence of 
masseless Dirac fermions with a Berry phase 7 = 0. Nev- 
ertheless, there is compelling evidence from far-infrared 
absorption, for the existence of Dirac fermions at the H 
point in graphite Far infrared measurements probe 
carriers in the very close vicinity of the H point where 
there is a maximum in the joint (initial and final) density 
of states. Transport measurements however, are sensitive 
to the density of states at Ep, which is modulated with 
increasing magnetic field, as the Landau bands cross the 
Fermi energy. For holes, maxima in the density of states 
correspond to Landau bands crossing Ep for k z < 0.5, 
away from the H point, where the dispersion is no longer 
linear and a priori there is no reason to expect the car- 
riers to behave as Dirac fermions. 

It is interesting to compare the data with the predic- 
tions of the SWM 3D band structure model with its seven 



> 

CD 

LU 




> 

CD 

E 

LU 



— i //'/ 
-(b) W/y 














. 1 


- . \A NT- 





4 6 
B(T) 



10 



FIG. 3: (color online) (a) Electron (increasing with B) and 
hole (decreasing with B) Landau bands (dE/dk z = 0) as a 
function of the magnetic field (calculated from SWM model 
for B > 0.4 T). The crossing of the electron and hole bands 
with the Fermi energy is in nearly perfect agreement with the 
measured magnetic field position of the electron (dotted lines) 
and hole (dashed lines) features, b) For B > 2 T, the Fermi 
energy is calculated self consistently, keeping the sum of the 
electron and hole concentrations constant. 



tight binding parameters 70, ••■,76- When the parameter 
73 is taken into account the magnetic field Hamiltonian 
has infinite order. This was numerically reduced to a 
600 x 600 matrix for the exact diagonalization procedure. 
A maximum in the conductivity is expected when there 
is a maximum in the density of states at the Fermi level. 
This can be found, for a given magnetic field, by looking 
for a Landau band with a slope dE / dk z = at Ep, which 
we refer to as a Landau band 'crossing' the Fermi level. 
Using the tight binding parameters of Ref. [55] as a start- 
ing point, this procedure was repeated until a satisfactory 
agreement with the data was obtained. The tight binding 
parameters found are given in Table [TJ While we are un- 
able to fit our data with exactly the same tight binding 
parameters as in Ref. |25j . the values we find are nev- 
ertheless not significantly different. Moreover, the pre- 
dicted (26] effective mass, m* — 4fi 2 7i/3ag7o = 0.054m e 
where eto = 0.246 nm is the in-plane lattice constant, cal- 
culated using our values for 70 and 71, is in good agree- 
ment with the accepted value [55] . 

For magnetic fields below 2 T it is a good approxi- 
mation to assume that the Fermi level is constant. The 
electron and hole Landau bands (solid lines), calculated 
using the parameters in Table [T] are plotted in Fig. [3^ a) 
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This work 


Ref. [25] 


7o (eV) 


3.37 ±0.02 


3.16 ±0.05 


7i (eV) 


0.363 ±0.05 


0.39 ±0.01 


72 (eV) 


-0.0243 ± 0.001 


-0.02 ±0.002 


73 (eV) 


0.31 ±0.05 


0.315 ±0.015 


74 (eV) 


0.07 ±0.01 


0.044 ± 0.024 


75 (eV) 


0.05 ±0.01 


0.038 ± 0.005 


76 = A (eV) 


-0.007 


-0.008 ±0.002 


E F (eV) 


-0.0287 


-0.024 ± 0.002 


9s 


2.4 ±0.1 




n (cm~ 3 ) 


-(2.4 ±0.4) x 10 17 





for massive fermions. No evidence for Dirac fermions 
at the H point is observed in the transport data. This 
can be understood, since transport is dominated by the 
crossing of the Landau bands at the Fermi level, corre- 
sponding to dE/dk z = 0, which occurs away from the 
H point (k z — 0.5), where the carriers are indeed Dirac 
fermions [HI2T]. 

This work has been partially supported by ANR con- 
tract PNANO-019-06. We acknowledge useful discus- 
sions with LA. Luk'yanchuk and Y. Kopelevich. 



TABLE I: Summary of the SWM tight binding parameters 
found in this work and compared to the values given in 
Ref. [25]. 



for magnetic fields below 2 T. The vertical broken lines 
indicate the observed electron (dotted) and hole (dashed) 
minima in AR XX (maxima in a xx ). The agreement be- 
tween the magnetic field position of the Landau bands 
crossing the Fermi level and the features in the transport 
data is remarkable. 

At higher magnetic fields, as graphite approaches the 
quantum limit, the Fermi energy is no longer constant 
as carriers are transferred between the electron and hole 
pockets. This is the reason for the considerable deviation 
from 1/B periodicity observed at high magnetic fields in 
Fig. (2jc). Nevertheless, as can be seen in Fig. [3jb) , the 
SWM model can correctly predict the magnetic field po- 
sition of the features provided the movement of Ep is 
taken into account. Here the Fermi level has been cal- 
culated self-consistently assuming the sum of the elec- 
tron and hole concentrations is constant, n — p = no. 
The electron concentration corresponds to the number 
of states in partially filled bands below the Fermi en- 
ergy, and the hole concentration to those above the 
Fermi energy. In order to fit the low-field data, we 
have used no = —2.4 x 10 17 cm~ 3 , assuming that un- 
der neutrality conditions, n = p = 8 x 10 18 cm~ 3 , with 
n = 64 7l | 72 |/(9^ 2 7o 2 a2 Co ) [25]. 

To reproduce the spin splitting in the high magnetic 
field data, a g-factor g s — 2.4 is required. In graphite, the 
g-factor cannot be reliably estimated from the separation 
AB of the spin split features since the Fermi energy is 
moving with field [7]. While this does not noticeably 
shift the orbital features below B = 2 T, the shift is 
significant compared to the spin gap. The value g s = 2.4 
should be considered as a lower limit. Any significant 
Landau level broadening, neglected in our model, would 
reduce the movement of the Fermi energy, and therefore 
increase the value of g s required to fit the data. 

To conclude, low temperature magnetotransport data 
of natural graphite and HOPG can be fully explained us- 
ing the Slonczewski-Weiss-McClure tight binding model 
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